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Abstract 

We present some Maple procedures using the GrTensorll package 
for teaching purposes in the study of the canonical version of the 
general relativity based on the ADM formalism 

1 Introduction 

The use of computer facilities cam be an important tool for teaching 
general relativity. We have experienced several packages of procedures, 
(in REDUCE + EXCALC for algebraic programming and in Math- 



ematica for graphic visualizations) which fulfill this purpose ([1C]). 
In this article we shall present some new procedures in MapleV using 
GrTensorll package ( |TT[| ) adapted for the canonical version of the gen- 
eral relativity (in the so called ADM formalism based on the 3+1 split 
of spacetime) . This formalism is widely used ( || , §§ ) in the last years 
as a major tool in numerical relativity for calculating violent processes 
as, for example the head-on collisions of black holes, massive stars or 
other astrophysical objects. Thus we used these computer procedures 
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in the process of teaching the canonical formalism as an introductory 
part of a series of lectures on numerical relativity for graduated stu- 
dents. The next section of the article presents shortly the notations 
and the main features of the canonical version of the general relativity. 
Early attemps in using computer algebra (in REDUCE) for the ADM 
formalism can be detected in the literature (H , || , jj]]) . Obviously we 
used these programs in producing our new procedures for Maple + 
GrTensorll package, but because there are many specific features we 
shall present in some detail these procedures in the section 3 of the 
article. The last section of the article is dedicated to the conclusions 
pointed out by running the Maple procedures presented here and some 
future prospectives on their usage toward the numerical realativity. 

2 Review of the canonical formalism 
of general relativity 

Here we shall use the specific notations for the ADM formalism Q , |@] ; 
for example latin indices will run from 1 to 3 and greek indices from 
to 3. The starting point of the canonical formulation of the general 
relativity is the (3+l)-dimensional split of the space-time produced 
by the split of the metric tensor : 



/ (4) 9oo (4) <^\ (N k N k -N 2 Nj 
V (4) <fc ^9ij ) \ N { 9ij 



where is the riemannian metric tensor of the three-dimensional 
spacelike hypersurfaces at t = const, which realize the spacetime 
foliation. Here N is the "lapse" function and N l are the components 
of the "shift" vector g. 

The Einstein vacuum field equations now are (denoting by " •" the 
time derivatives) : 

gij = 2Ng- 1 / 2 [ nj - \gi^ k k ] + N ifj + N 3/i (2) 
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+ [^N m ] /m - N* /m ^ - Ni /m7 r mi (3) 

where it 13 are the components of the momenta canonically conjugate 
to the gij's. 

In the above formulas we denoted by "/" the three-dimensional 
covariant derivative defined with g^ using the components of the three- 
dimensional connection B : 

= -j9 m {9mi,h + 9mk,j - 9jk,m) (4) 
The Ricci tensor components are given by 

^H? ~~ 1 ij,k 1 ik,j ~r 1 ij 1 km 1 im 1 jk \°) 

The initial data on the t = const, hypersurface are not independent 
because they must satisfy the constraint equations, which complete 
the Einstein equations 

n = -Vg{R + g-^l^k? - = o (6) 

H l = -2TT lj /j = (7) 

where Ti is the super-hamiltonian, H % the super-momentum and g is 
the determinant of the three-dimensional metric tensor g^ . 

The action functional in Hamiltonian form for a vacuum space 
-time can thus be written as ( , JH ) : 

S = J dtj (ir ij gij - NH — NiH^uWu 3 (8) 

where the oj 1, s are the basis one- forms. Thus the dynamic equations 
(^) and (||) are obtained by differentiating S with respect to the canon- 
ical conjugate pair of variables (vr*-', gkm)- 



3 Maple + GrTensorll procedures 

Here we shall describe briefly the structure and the main features 
of the Maple procedures for the canonical formalism of the general 
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relativity as described in the previous section. Two major parts of 
the programs can be detected : one before introducing the metric of 
the spacetime used (consisting in several definitions of tensor objects 
which are common to all spacetimes) and the second one, having line- 
commands specific to each version. 

The first part of the program starts after initalisation of the GrTen- 
sorll package (grtw();) and has mainly the next lines : 

> grdef (<tr := pi{~i i}'); 

> grdef ( 'haO : =-sqrt (detg) * (Ricciscalar+ 

(l/detg)*((l/2)*(tr)~2-pi{i j}*pi{ "i ~j }))'); 

> grdef (<ha{ ~i }:=-2*(pi{ "i ~j ;j>-pi{ "i ~j }*Chr{ p j "p })'); 

> grdef (<derge{ i j }: =2*N(x,t) * (detg) ~ (-1/2) *(pi{ i j } - 

(l/2)*g{ i j}*tr)+Ni{ i ;j > + Ni{ j ;i }<); 

> grdef ('Ndd{ "m j > : = Nd{ "m ; j > ' ) ; 

> grdef ( 'bum{ "i "j "m}:=pi{ "i "j }*Ni{ "m }'); 

> grdef ('bla{ "i "j }:=bum{ "i "j "m ;m }'); 

> grdef ( 'derpi{ "i "j }:= 

-N(x,t)*(detg) " (1/2)*(R{ "i "j }- (l/2)*g{ "i "j >*Ricciscalar)+ 
(l/2)*N(x,t)*(detg)-(-l/2)*g{ ~i ~j }*(pi{ "k "1 }*pi{ k 1 }- 
(l/2)*(tr)-2)-2*N(x,t)*(detg)"(-l/2)*(pi{ ~i "m >*pi{ "j m >- 
(l/2)*pi{ "i "j }*tr)+ (detg)" (1/2) *(Ndd{ "i "j >-g{ "i "j }* 

Ndd{ "m m >) + bla{ "i "j > - Ni{ ~i ;m >*pi{ "m "j }- 

Ni{ "j ;m }*pi{ ~m "i >') ; 

Here haO and ha{ i } represents the superhamiltonian and the 
supermomentum as defined in eqs. @ and (0) respectively and tr 
is the trace of momentum tensor density ir 11 - which will be defined 
in the next lines of the program. Here N(x,t)represents the lapse 
function N. Also, derge{ i j } represents the time derivatives of the 
components of the metric tensor, as defined in eq. (|2|) and derpi{ * i 
j } the time derivatives of the components of the momentum tensor 
■k'^ as defined in eq. (||). 

The next line of the program is a specific GrTensorll command 
for loading the spacetime metric. Here Maple loads a file (previously 
generated) for introducing the components of the metric tensor as 
functions of the coordinates. We also reproduced here the output of 
the Maple session showing the metric structure of the spacetime we 
introduced. 

> qload('Cyl_din') ; 
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Default spacetime = Cyl_din 
For the Cyl_din spacetime: 
Coordinates 
x(up) 

a 

x = [x, y, z] 
Line element 
2 2 
ds = exp(gamma(x, t) - psi(x, t)) d x 

2 2 2 

+ R(x, t) exp(-psi(x, t)) d y + exp(psi(x, t)) d z 

As is obvious we introduced above the metric for a spacetime with 
cylindrical symmetry, an example we used for teaching purposes being 
a well known example in the literature (H). In natural output this 
metric has the form : 



9ij 



( ei~^ 

R 2 e~^ j (9) 
v 



in cylindrical coordinates x,y,z with x S [0, oo), y G [0, 2tt), z G 
(— oo, +oo) where R, if) and 7 are functions of x and t only. 

After the metric of the spacetime is established the next sequence 
of the programm just introduce the components of the momentum 
tensor tt ij as 

> grdef ('Nd{ " m > := [dif f (N(x,t) ,x) , 0, 0]'); 

> grdef('Ni{ "i > := [Nl(x,t), M2(x,t) , N3(x,t)]'); 

> grdef ('vil{~i}:= [pig(x, t) *exp(psi (x,t) -gamma (x,t) ) ,0,0]'); 

> grdef ('vi3{~i> 

:=[0,0,exp(-psi(x,t))*(pig(x,t)+(l/2)*R(x,t)*pir(x,t)+ 

pip(x,t))] ') ; 

> grdef ( ( vi2{"i>:=[0, (2*R(x,t))"(-l)*pir(x,t)*exp(psi(x,t)) ,0] ') ; 

> grdef Cpi{ "i ~j > : = 

vil{ "i >*kdelta{"j $x}+vi2{ ~i }*kdelta{ "j$y }+ 
vi3{ "i }*kdelta{~j $z}'); 

> grcalc(pi(up,up) ) ; 

> grdisplay(pi(up,up) ) ; 
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Here Ni{ " i } represents the shift vector N' 1 and the other objects 
(Nd, vil, vi2 and vi3) represent intermediate vectors defined in order 
to introduce the momenum pi{ " i " j } having the form : 



7f 




7r 7 e 








,1/1-7 




(10) 



In the program we denoted 7r 7 , ttr and tt^ with pig, pir and pip, 
respectively. The momentum components are introduced in order that 
the dynamic part of the action of the theory be in canonical form, that 
is : gijiT 13 = 7r 7 7 + -K^pip + ttrR. The next lines of the programm check 



if this condition is fullfiled : 

> grdef ('del{ i > : = [dif f (grcomponent (g(dn,dn) , [x,x] ) ,t) , 0,0] ' ) ; 

> grdef ('de2{ i >:= [0 , dif f (grcomponent (g(dn,dn) , [y,y] ) ,t) ,0] ') ; 

> grdef ( c de3{ i }:= [0 , ,diff (grcomponent (g(dn,dn) , [z,z] ) ,t)] ') ; 

> grdef (<ddgt({ i j }:= 

del{ i }*kdelta{j $x>+de2{ i }*kdelta{ j$y }+ 
de3{ i }*kdelta{ j $z}'); 

> grcalc (ddgt (dn , dn) ) ; 

> grdef ('act: =pi{ ~i "j >*ddgt{ i j }'); 

> grcalc (act) ; gralter (act, simplify) ; grdisplay (act) ; 

By inspecting this last output from the Maple worksheet, the user can 
decide if it is necessary to redifine the components of the momentum 
tensor or to go further. Here the components of the momentum tensor 
were calculated by hand but, of course a more experienced user can try 
to introduce here a sequence of commands for automatic calculation of 
the momentum tensor components using the above condition, through 
an intensive use of solve Maple command. 

Now comes the must important part of the routine, dedicated to 
calculations of different objects previously defined : 

> grcalc (haO) ; gralter (haO, simplif y) ; 

> grdisplay (haO) ; 

> grcalc (ha (up) ) ; gralter (ha(up) , simplif y) ; 

> grdisplay (ha(up) ) ; 

> grcalc (derge (dn, dn) ) ; gralter (derge (dn, dn) , simplify) ; 

> grdisplay (derge (dn, dn) ) ; 

> dl : =exp(-psi (x, t) ) *grcomponent (derge (dn,dn) , [z ,z] )+exp(psi(x, t) - 
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gamma (x,t))*grcomponent(derge(dn,dn) , [x,x] ) ; 

> simplify(dl) ; 

> d2:=(l/(2*R(x,t)))*exp(psi(x,t))*grcomponent(derge(dn,dn) , [y,y] )+ 

(l/2)*R(x,t)*exp(-psi(x,t))*grcomponent(derge(dn,dn) , [z,z] ) ; 

> simplify(d2) ; 

> d3:=exp(-psi(x,t))*grcomponent(derge(dn,dn) , [z,z] ) ; 

> simplify (d3) ; 

> grcalc (derpi (up ,up) ) ; gralter(derpi(up, up) , simplify) ; 

> grdisplay (derpi (up ,up) ) ; 

> fl := exp(gamma(x,t)-psi(x,t))*grcomponent(derpi(up,up) , [x,x] )- 

pig(x,t)*(d3-dl) ; 

> simplif y(f 1) ; 

> f2:= 2*R(x,t)*exp(-psi(x,t))*grcomponent(derpi(up,up) , [y,y] )+ 

( 1/R (x , t ) ) *d2*pir (x , t ) -pir (x , t ) *d3 ; 

> simplify(f2) ; 

> f3 := exp (psi(x,t))*grcomponent (derpi (up, up) , [z,z] )+d3*(pig(x,t)+ 

(l/2)*R(x,t)*pir(x,t)+pip(x,t))-f l-(l/2)*R(x,t)*f2- 

(l/2)*pir(x,t)*d2; 

> simplify(f3) ; 

This is a simple series of alternation ofgrcalc, gralter and grdisplay 

commands for obtainig the superhamiltonian, supermomentum and 
the dynamic equations for the theory, dl ... d3 and fl ... f3 are 
the time derivatives of the dynamic variables, 7, R, ip and tt 7 , ttr, 
■fr^ respectively. Denoting with "/" the derivatives with respect to r 
we display here the results for the example used above (cylindrical 
gravitational waves) : 

H° = e^(2R" - RW + \w?R - tt.ttr + ^M 2 ) = 

H 1 =H r = e^~ 1 (-2it' 1 + 7'7r 7 + R'ir R + ip'n.p) = ; H 2 = H 3 = 

7 = N 1 -/' + 2N 1 ' - e^NiTR ; R = N 1 R' - Nit^ 

R 

^ = N 1 ^ + JV 1 '^ - e^iR'N' + -R'tp'N - ^' 2 RN + -N^ttr - ^MrJ) 



7 




These are the well-known results reported in (Q) or (||). 

One of the important goals of the canonical formalism of the gen- 
eral relativity (which constitutes the "kernel" of the ADM formalism) 
is the reductional formalism. Here we obtain the true dynamical sta- 
tus of the theory, by reducing the number of the variables through 
solving the constraint equations. This formalism is applicable only 
to a restricted number of space-time models, one of them being the 
above cylindrical gravitational waves model. Unfortunately only a 
specific strategy can be used in every model. Thus the next lines of 
our program must be rewritten specifically in every case. Here, for 
teaching purposes we present our example of cylindrical gravitational 
wave space-time model. Of course we enccourage the student to apply 
his own strategy for other examples he dares to calculate. 

In our example of cylindrical gravitational waves, the reductional 
strategy as described in (H) starts with the usual rescaling of TL and 



W to H and W by 

H = e^H ; N = eVjV ; H 1 = e^H 1 ; N 1 = e^N 



wich produce the next sequence of Maple+GrTensorll commands : 

> grdef ( ' ahaO : =sqrt ( exp (gamma (x, t) -psi(x,t) ) ) *haO' ) ; 

> grdef ('aha{ "j } := exp (gamma (x,t)-psi(x,t))*ha{ "j }'); 

> grdef ( ' an: =sqrt (exp(psi (x,t) -gamma (x, t) ) ) *n(x, t) ' ) ; 

> grdef ('ani{ "i } := exp(psi(x,t)-gamma(x,t))*ni{ "i }'); 

The canonical transformation to the new variables, including Kuchar's 
"extrinsic time", defined by : 




oo 




R = R 



tt r + [In ( 



R' + T 1 
R! -V 



)]' 



are introduced with : 
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> pig(x,t) :=-diff (T(x,t) ,x) ; 



> pir(x,t) :=piR(x,t) - diff (ln((diff (R(x,t) ,x)+diff (T(x,t) ,x))/ 

(diff (R(x,t) ,x)-diff (T(x,t) ,x))) ,x) ; 

and specific substitutions in the dynamic objects of the theory : 

> grmap(haO, subs , dif f (gamma(x,t) ,x) = 

diff( ln( (diff (R(x,t) ,x))"2- (dif f (T(x,t) ,x) ) ~2 ) ,x)-piT(x,t) , 'x' ) ; 

> grcalc(haO) ; gralter (haO, simplify) ; 

> grdi splay (haO) ; 

> grmap(ha(up) , subs , dif f (gamma(x,t) ,x)=dif f ( ln( (dif f (R(x,t) ,x) ) "2- 

(diff (T(x,t) ,x))~2 ) ,x)-piT(x,t) , 'x') ; 

> gralter (ha(up) , simplify) ; 

> grdisplay (ha(up) ) ; 

> grcalc (ahaO) ; 

> grmap(ahaO, subs , dif f (gamma(x,t) ,x)=dif f ( ln( (dif f (R(x,t) ,x) ) "1- 

(diff (T(x,t) ,x))-2 ) ,x)-piT(x,t) ,'x'); 

> gralter (ahaO, simplif y, sqrt) ; 

> grdisplay (ahaO) ; 

> grcalc (aha(up) ) ; 

> grmap(aha(up) , subs , dif f (gamma(x,t) ,x)=dif f ( ln( (dif f (R(x,t) ,x) ) "2- 

(diff (T(x,t) ,x))-2 ) ,x)-piT(x,t) ,'x'); 

> gralter (aha(up) , simplify) ; 

> grdisplay (aha(up) ) ; 

> grmap(act, subs , dif f (gamma(x,t) ,x)=dif f ( ln( (dif f (R(x,t) ,x) ) "2- 

(diff (T(x,t) ,x))"2 ) ,x)-piT(x,t) ,'x'); 

> grcalc(act); grdisplay (act) ; 

Thus the action yields (modulo divergences) : 

/oo rco _ _ _ _ 

dt / dr(Il T f + U R R + 7T^V> + tt x x — NH — N{H l ) 
-oo J 

where : 

H = R'U T + T'U R + X -R- x ^% + X -Ri)' 2 + \R-^ 2 X + R X ' 2 
H 1 = T'U T + R'U R + 4)'^ + x'vr x 
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Solving the constraint equations H = and H = for Hy and 11^ 
and imposing the coordinate conditions T = t and R = r we obtain 
finally : 

/+oo r+oo ^ 

dT / dR[ir^, T + 7t xX ,t - AR~^l + ty% + R*l + R^X' 2 )} 
-oo JO ^ 

from the next sequence of programm lines : 

> R(x,t):=x; T(x,t):=t; grdisplay (ahaO) ; 

> solve (grcomponent (ahaO) ,piT(x,t)) ; 

> piT(x,t):= -l/2*(x~2*diff (psi(x.t) ,x) ~2+pip(x,t) ~2)/x; 

> eval(piR(x,t)) ; 

> piR(x,t) :=-diff (psi(x,t) ,x)*pip(x,t) ; piR(x,t) ; 

> grdisplay (ahaO) ; grdisplay (aha(up) ) ; 

> piT(x,t) ; 

2 /d \2 2 
x | — psi(x, t) | + pip(x, t) 
\dx / 
- 1/2 

x 

> piR(x,t) ; 

/d \ 
-| — psi(x, t) | pip(x, t) 
\dx / 

> grcalc(act); grdisplay (act) ; 

For the Cyl_din spacetime: 

act 

/d \ 
act = | — psi(x, t) | pip(x, t) 
\dt / 

> grdef ('Action:=act+piT(x,t)*diff (T(x,t) ,t)+piR(x,t)*dif f (R(x,t) ,t) ') ; 

> grcalc (Action) ;gralter(Action, factor, normal, sort, expand) ; 
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> grdisplay (Action) ; 

For the Cyl_din spacetime: 



Action 



/d \2 /d \ 

Action = - 1/2 x | — psi(x, t) | + | — psi(x, t) | pip(x, t) 
\dx / \dt / 

2 

pip(x, t) 
- 1/2 

x 

> grdef ('Ham:=piT(x,t)*diff (T(x,t) ,t)+piR(x,t)*diff (R(x,t) ,t) ') ; 

> grcalc(Ham) ; gralter (Ham, expand) ; 

> grdisplay (Ham) ; 

For the Cyl_din spacetime: 
Ham 

2 

/d \2 pip(x, t) 

Ham = - 1/2 x |— psi(x, t) | - 1/2 

\dx / x 



4 Conclusions. Further improuvements 

We used the programms presented above in the computer room with 
the students from the graduate course on Numerical Relativity. The 
main purpose was to introduce faster the elements of the canonical 
version of relativity with the declared objective to skip the long and 
not very straitforward hand calculations necessary to process an entire 
example of spacetime model. We encouraged the students to try to 
modify the procedures in order to compute new examples. 

The major conclusion is that this method is indeed usefull for an 
attractive and fast teaching of the methods involved in the ADM for- 
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malism. On the other hand we can use and modify these programs for 
obtaining the equations necessary for the numerical relativity. In fact 
we intend to expand our Maple worksheets for the case of axisymmet- 
ric model (used in the numerical treatement of the head-on collision 
of black- holes). Of course, for numerical solving of the dynamic equa- 
tions obtained here we need more improuvements of the codes for 
paralel computing and more sophisticated numerical methods. But 
this will be the object of another series of articles. 
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